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Abstract 

We present a new phenomenological analysis of the inclusive rare decay B — > XJ + l~ . 
In particular, we present the first calculation of the NNLL contributions due to the 
leading two-loop matrix elements, evaluated for arbitrary dilepton invariant mass. This 
allows to obtain the first NNLL estimates of the dilepton mass spectrum and the lepton 
forward-backward asymmetry in the high Mf+£- region, and to provide an independent 
check of previously published results in the low M| + ^_ region. The numerical impact of 
these NNLL corrections in the high-mass region {Mj +i - > 14.4 GeV 2 ) amounts to -13% 
in the integrated rate, and leads to a reduction of the scale uncertainty to ±3%. The 
impact of non-perturbative contributions in this region is also discussed in detail. 
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1 Introduction 



It has been more than a quarter of a century since the basic blocks were laid to build up the 
Standard Model (SM). Since then a wealth of high precision calculations and measurements has 
been performed and compared. The SM passed all the tests. And yet, it is universally accepted 
that there must be something more than what has been so brilliantly conceived and verified. Our 
undaunted quests to understand fundamental problems such as the flavour families, the dark matter 
or quantum gravity, have led to several proposals of new models. Short of direct evidences about 
the new degrees of freedom of the theory, it remains a fact that some of the most constructive 
and definitive work to search for extensions of the SM lies in the cold hard world of higher and 
higher precision tests: the systematic search for discrepancies between theoretical estimates and 
experimental data in quantities particularly sensitive to the symmetry structure of the model. 

Various possibilities of precision tests have been considered, some of which derive from flavour- 
changing neutral-current processes (FCNC). The key argument here is that neutral flavour- violating 
processes can occur only via loops, within the SM, and leads to tight constraints on possible new 
sources of flavour-symmetry breaking terms, even if these appear well above the electroweak scale 
(see e.g. Refs. 0121 El)- The new B factories, BELLE and BABAR, are providing us now with more 
and more high statistics data, which will carry FCNC tests to the next precision level. For our part, 
we shall now focus on an inclusive process of the type B — > X s £ + £~ , where in principle X s stands for 
anything with a single strange quark. Since the emphasis is on precision, we must work in a region 
where uncertainties can be minimized or at least reliably bounded. As long as the (sub)energy 
scales are much larger than the characteristic Aqcd we can treat the problem on hand primarily 
at a parton level, where perturbative calculations are justifiably deployed. Non-perturbative effects 
either must be under control or even can be added with confidence. These requirements call for our 
avoiding the regions where the dilepton invariant mass is near the cc resonances. 

After such considerations, it is natural to divide the B — > X s £ + £~ data into two distinct sets, 
depending on the squared invariant mass (Mj +i - = q 2 ) of the dilepton pair: 

1 GeV 2 < q 2 < 6 GeV 2 [low]; q 2 > 14.4 GeV 2 [high]. 

These two windows are in fact complementary to each other, both experimentally and theoretically, 
on the one hand because of identification and detection efficiency and event rates, on the other 
because of non-perturbative and perturbative effects, as well as parametric uncertainties. We shall 
give results that will cover these two regions. 

It was recognized some time ago that, because of the mixing structure and because of the very 
heavy top mass, rates for processes such as B — > X s ^ and B — > X s £ + £~ could be sizable. It was 
also realized that these processes are ideal candidates for a short- distance analysis based on the 
operator product expansion. Partonic QCD corrections to the pure electroweak amplitudes can 
be investigated systematically, via renormalization group improvement, and are found to change 
the rates significantly. With the help of the heavy-quark expansion one can also reliably assess 
non-perturbative contributions, which are found to be small for sufficiently inclusive observables. 
As one may infer from above, there are essentially three ingredients which go into a high-accuracy 
calculation of these FCNC processes. The first one is the perturbative evaluation of the partonic 
amplitudes at the electroweak scale, which can be translated into the initial conditions for an 
effective theory description based on a local Hamiltonian, TC e s = Ci Oi, of dimension-6 operators. 
The second ingredient deals with the subsequent evolution of TC e s down to the actual scale of 
the process, namely the running of the effective coupling constants Cj - which tentamounts to 
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dressing them with QCD - from Mw,t to mj via renormalization group improvement. The last 
step concerns the evaluation of a combination of process-dependent matrix elements corresponding 
to the specific observable: a combination of the type / dT\{X s £ + £~\ J2C{Oi\B)\ 2 in our case. This 
last step includes hard QCD corrections, which again can be computed perturbatively, and the 
attendant bremsstrahlung and soft gluon effects, since QCD is a massless gauge theory. It also 
includes the non-perturbative hadronization of b and s quarks, which can be analyzed with the 
help of the heavy-quark expansion. Each of the three steps must be taken to matching orders of 
accuracy in powers of the strong coupling constants a s , with renormalization group resummations 
whenever necessary. Within the present paper, our central focus will be on the last step (hadronic 
matrix elements) and we will give expressions for the differential rate in the dilepton mass and 
the integrated rates over the two windows. We already published results for the forward-backward 
asymmetry [3J |3] , where we detailed our calculation of the bremsstrahlung contributions. 

One may view B — ► X s £ + £~ as an extension of B — > X^, in which the virtual photon can 
take on a whole spectrum of mass. This picture is somewhat simplistic and rather misleading. 
Indeed, to describe the former process we need extra electroweak operators, which are sensitive 
to different aspects of the mechanisms of electroweak- and flavor-symmetry breaking. This is a 
very important feature, which could allow to detect a non-standard effect in B — > X s £ + £~ even in 
absence of deviations from the SM in B — > X s ^y. On the other hand, there are certainly several 
common points in the theoretical analysis of these two processes. In particular, similarly to the 
B — > Xgj case, also in B — > X s £ + £~ sizable corrections to the pure electroweak amplitude are 
induced by the leading four-quark operators Oi,2- Thus the matrix elements (s£ + £~\Oi,2\b) must 
be evaluated accurately, beyond the first non-trivial order. The corresponding Feynman diagrams 
are one-gluon-corrections to the basic one-loop diagrams with a charm quark contracted. These 
correspond to two-loop integrals with three relevant scales: m&, m c and q 2 (not simply two mass 
scales as in B —>■ Xgj). So far these diagrams have been analyzed only in Ref. |BJ, by means of a 
mass and momentum double expansion method. By its very nature of expansion, the method of 
Ref. [U] is not applicable when the dilepton mass approach the cc threshold (q 2 ~ 4m 2 ) or becomes 
larger. We have applied a semi-numerical method, which first converts analytically all two-loop 
integrals into a standard set of integrals and then performs a rapid numerical integration over a 
set of Feynman parameters. This method is very accurate and works over any physical kinematical 
range, thus we are able to predict the partial decay rate not only in the low-mass window (as in 
Ref. [6]), but also above the cc threshold. As a result, our present work provides both a detailed 
and accurate check of the results of Asatrian et al. in the low mass window and, at the same time, 
allow us to present the first NNLL results for the high-mass window. 

The NNLL perturbative calculations will bring predictions of higher precision only if we can have 
the same level of confidence when estimating the non-perturbative effects. These are divided into 
hadronization of the external b and s quarks, and the residual long-distance effects due to the tails 
of the narrow cc resonances. In both cases the corrections can be analyzed using appropriate heavy- 
mass expansions and the resulting uncertainties will be shown to be under reasonable control in 
both windows. In the high-mass range non-perturbative uncertainties turns out to be the dominant 
source of theoretical errors. However, as we shall show, a considerable reduction of this uncertainty 
can be expected in the near future with a better knowledge of universal non-perturbative parameters 
from other processes. 

The experimental situation regarding the inclusive decay B — > X s £ + £~ is as follows: BELLE and 
BABAR have already obtained clear evidences (~ 5a) of this transition, quoting two measurements 
which are compatible with each other and with the SM expectation |S] . Both results are based 
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on a semi-inclusive analysis: the hadronic system X s is reconstructed from a kaon plus to 4 pions 
(at most one 7r°). The signal characteristics is determined by modeling the invariant mass Mx s 
spectrum using the phenomenological model first proposed in [9 . The reconstruction efficiencies of 
the signal are determined by the MC samples based on this model, which leads to a large fraction 
of the present systematic uncertainty. The overall uncertainty of these first measurements of the 
inclusive decay rate is still at the 30% level, but substantial improvements can be expect in the near 
future. 

The plan of this paper is as follows: in the next section we shall present a more thorough 
discussion of the effective-theory approach to inclusive FCNC b decays, briefly reviewing what is 
already known in the literature. In Section 3, we shall describe in more detail the method we have 
used to calculate the two- loop matrix elements (s£ + £~\Oi t 2\b). We shall give plots of these matrix 
elements and compare them with those obtained by Asatrian et al.: besides a very good agreement 
within the low-g 2 regime, they also show how and when the expansion method fails as we approach 
the cc threshold; our new results on the necessary one-loop counterterms are also presented. In 
Section 4 we present our new results for the matrix element of 8 over the whole dilepton invariant 
mass spectrum. Section 5 is reserved for a concise analysis of the non-perturbative effects which 
also includes 1/mf contributions. .In Section 6 we finally present our phenomenological analysis: we 
quote results for the integrated rates over the two disjoint dilepton mass spectrum windows, for the 
integrated lepton forward-backward asymmetry in the high mass region and for the position zero 
of the asymmetry. A proper assignment of the theoretical uncertainties associated to these results 
will be presented. A few appendices have been prepared to facilitate the reading of the main text. 



2 Theoretical framework 
2.1 Effective field theory 

Within inclusive B decays, such as B — > A s 7 and B — > X s £ + £~, short- distance QCD corrections 
are sizable and comparable to the pure electroweak contributions. They stem from the evolution 
of the system from a large scale M^ eaxj = O(Mw), where the weak interaction acts, to the decay 
energy ~ m&, resulting in large logarithms of the form a™(m&) log m (mb/Mh ea vy), where m < n 
(with n = 0, 1, 2, ...). The most suitable framework for their necessary resummations is an effective 
low-energy theory with five quarks, obtained by integrating out the heavy degrees of freedom. 
Renormalization-group (RG) techniques are used to organize the resummation of the series in leading 
logarithms (LL), next-to- leading logarithms (NLL), and so on: 

<(m 6 ) \og n (m b /M) [LL], a n s + \m b ) log" (m b /M) [NLL],... (2.1) 

The effective five-quark low-energy Hamiltonian relevant to the partonic process b — > s£ + £~ can be 
written as 

H eff = ^V t *V tb £ GO*, M heavy ) Ofo) , (2.2) 
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where 



Oi = (s lfl T a P L c){c^T a P L b), 2 
5 = (s ltlll/lp P L b)E q (qri u Yq), O, 

lD7T / 
lD7T z 

define the complete set of relevant dimension-6 operators; Cj(/x, Mh eavy ) are the corresponding Wil- 
son coefficients. As the heavy fields are integrated out, the top-, W-, and Z-mass dependence is 
contained in the initial conditions of these Wilson coefficients, determined by a matching procedure 
between the full and the effective theory at the high scale (Step 1). By means of RG equations, the 
Cj(/x, M heavy ) are then evolved to the low scale (Step 2). Finally, the QCD corrections to the matrix 
elements of the operators are evaluated at the low scale (Step 3). 

Compared with the effective Hamiltonian relevant to b — > 57, Eq. (J2.3|) contains additional 
operators Og and O 10 which are of order a em . The first large logarithm of the form \og(m b /M w ) 
already appears without gluons, because the operator C 2 mixes into 0$ at one loop. It is then 
convenient to redefine the magnetic, chromomagnetic and lepton-pair operators as follows [HH lllj: 

16tt 2 ~ q 2 ~ 

O i = —O i , C^jf^Q, (z = 7,...,10). (2.4) 

With this redefinition, one can follow the three calculational steps discussed above ITU} lllj. as in 
b — > 57. In particular, after the reshufflings in ()2.4|) the one-loop mixing of the operator O2 with 
O9 appears formally at order a s . At NNLL precision, one should in principle take into account the 
QCD two-loop corrections to (O9), the QCD one-loop corrections to (Or) and (O10), and the QCD 
corrections to the electroweak one-loop matrix elements of the four-quark operators. 



= (s 7m Plc) (c^P L b) , 

= (s^PL^Egiq^Taq), 

= (s lllll/lp T a P L b)Z q (qrYYT a q) , 

(2.3) 

= ^m b Oi)(Sa^T a P R b)G% i 
= (s^P L b) {lYld) 



2.2 Present status of the calculation 

Regarding the present status of these perturbative contributions to decay rate and FB asymmetry 
of B — > X s £ + £~ (for a recent review see j3]), we note that the complete NLL contributions to the 
decay amplitude can be found in [TUl HI] . Since the LL contribution to the rate turns out to be 
numerically rather small, NLL terms should be regarded an 0(1) correction to this observable. On 
the other hand, since a non-vanishing FB asymmetry is generated by the interference of vector 
(~ 07,9 ) and axial-vector (~ O10) leptonic currents, the LL amplitude leads to a vanishing result 
and therefore NLL terms become by default the lowest non-trivial contribution to this observable. 

For these reasons, a computation of NNLL terms in B — > X s £ + £~ is needed if we are to aim for 
the same numerical accuracy as achieved by the NLL analysis of B — > X^. Before we proceed any 
further, we would like to acknowledge the accumulated efforts by many groups, whose contributions 
greatly lessen the toil. Indeed, a large body of results for b — > 37 can be taken over and used in the 
NNLL calculation of B —* X s £ + £~ ^21 EE El E3 E2] • The necessary additional calculations, includ- 
ing the two-loop corrections presented in this paper, have been cross-checked and our contribution 
here is a part of that finalization process. 
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To begin, the full computation of initial conditions to NNLL precision was presented in Ref . jT7] . 
It removes the large matching scale uncertainty, which amounts to around 16% in the NLL approx- 
imation. 

Part of the additional three-loop mixings (Step 2), which were not known from the B — > X s ^ 
case, and which connect one of the four-quark operators C?i_6 to O7 and Og, have been reported 
recently in [18] : their effect leads to a 2% correction of the rate. The NNLL intra- mixings of the 
four-quark operators are still missing but are expected to yield an even smaller contribution. Since 
the FB asymmetry does not receive contributions from the term proportional to |(Og)| 2 , these 
mixings terms are not needed for a NNLL analysis of this observable. 

In Step 3 the most important NNLL contributions come from the two-loop matrix elements of 
the four-quark operators 0\ and Oi- They were calculated in [Hj, using Mellin-Barnes techniques 
similar to the ones originally used in the corresponding B — > X s 7 calculation ^Hj- These lead to 
a double expansion in m c /m b and q 2 /ml, where q 2 is the dilepton mass squared. Thus, the results 
in [Hj are only valid below the cc threshold within the dilepton mass spectrum. The calculation of 
the NNLL two-loop matrix elements advocated by us in this paper is based on a semi-numerical 
method. Since its validity spans over the whole dilepton mass spectrum, our work not only gives an 
independent check of the calculation in |Hj within the low q 2 window, but stands on its own merits 
to give results above the cc threshold. A complete NNLL calculation also requires the one-loop 
matrix element of the operator Og- This was done for small dilepton mass in We now present 
results that are valid for arbitrary dilepton mass, to be reported below. In 6 3 it was shown that 
within the integrated low-dilepton mass spectrum these NNLL matrix element contributions reduce 
the perturbative uncertainty (due the low-scale dependence) from ±13% down to ±6.5% and also 
the central value is changed significantly, ~ 14%. We shall report our findings for both windows 
in due course. The NNLL bremsstrahlung contributions and the corresponding IR virtual one-loop 
corrections have also been calculated for the dilepton mass spectrum (symmetric part) in |U HH| 
and for the FB asymmetry in [311201 EI] Their impact is significant; the zero of the FB asymmetry 
for example gets shifted by 11% by these contributions. 

For completeness, we note that within Step 3 a few pieces are still missing but their effects are 
estimated not to exceed 2%. In particular, a complete NNLL calculation of the B — > X s £ + £~ rate 
would require also the evaluation of the two-loop matrix element of the operator Og: its influence 
on the dilepton mass spectrum is expected to be very small, because it gets multiplied by a small 
leading Wilson coefficient Cg°\ In any event, similar to the missing entries of the anomalous- 
dimension matrix, this (scale-independent) contribution does not enter into the FB asymmetry at 
NNLL accuracy. The list of missing contributions includes also a calculation of two-loop matrix 
elements of the four-quark (penguin) operators 03,4,5,6- However, their net NNLL effect is most likely 
to be strongly suppressed by their small Wilson coefficients. This expectation is substantiated by 
the corresponding contributions to the B — > X s 7 decay branching ratio. The latter are shown to 
have an effect below 1% [22 J. 

In conclusion, the QCD NNLL calculation of the FB asymmetry is now fully complete, while 
the one for the dilepton mass spectrum distribution is on the verge of being completed. Our present 
work is a part of that endeavor. All other missing pieces can be estimated to be smaller than 2%. 
At this level of accuracy, other subleading effects may turn out to be more important. In particular, 
the uncertainties regarding input parameters should deserve a lot of attention, and further studies 
regarding higher-order electromagnetic effects are necessary. 
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2.3 Basic expressions 

In this subsection, we want to recapitulate some expressions and definitions of the basic observables. 
We normalize all by the semileptonic decay width in order to reduce the uncertainties due to bottom 
quark mass and CKM angles: 



F[b — > X c ev e 



(2.1) 



Here z = m 2 /m 2 (m& )C denote pole quark masses), f(z) = 1 — 8z + 8z 3 — z 4 — 12z 2 lnz is the 



phase-space factor and 



k{z) 



2a s (m b ) h(z) 



(2.2) 



3tT f(z) 

takes into account QCD corrections. The function h(z) has been given analytically in and is 
explicitly displayed in one of our appendices. The normalized dilepton invariant mass spectrum is 
then defined as 



R(s) 



(2.3) 



F(B -> X c ev) ' 

where s = (pi+ + pi-) 2 /m 2 . The other important observable is the forward-backward lepton asym- 
metry: 

A FB ( S ) = — — — / dcos9 e \ s ^sgn(cos^), (2.4) 

I [B — > X c ev) J -i as a cos 6^ 

where 0£ is the angle between the £ + and B momenta in the dilepton center-of-mass frame. It was 
shown in [21] that Afb{s) is equivalent to the energy asymmetry introduced in [2*5] . 

We also present here some useful formulae that will allow us to systematically take into account 
all corrections to these two observables at the partonic level beyond the NLL approximation: 



R(s) 



a 



2 T/*T/ 2 

V tb V ts 



47T 5 



V, 



cb 



1-8)2 r4 ( 1+ ?" c 



f{z)n(z 
+(l + 2s) 



r(s)\ 2 



at. 



71 



|c 9 ncw ( s )| 2 + \a 



™(s)\ 2 



a 



1 H r 99 (s) 



71 



+ 12 Re [Cr(s)Cr(sT] (l + ^r 79 ( S )) + 



A FB (s) 



3 «L 




2 (i--) 2 f 


4vr 2 




/(*)*(*) i 



S Re [C^ ew ( S )*C 9 new ( S )] 1 + -^r 910 (s 



7T 



+2Re[C^(s)*Cr w (s)} (l + ^tyio(s)) + ^fb(s)} 
In these expressions, we have introduced a new set of effective coefficients, defined as 

< 0) FP{ 8 ) + Ci 0) F?\s) + Cf (0) FP 
\c{ 0) Fl 9 \s) + cf F 2 (9) ( S ) + &f®Fi% 

c? cw (*) 



C? ew (s) 



i + -<7 7 ( s ) ) cf - ^ crFi'\s) + cr^» + cr u; *r(*) 

7T / 4 7T L 

l + —a 9 (s))Cf(s)- — 

71 J 4 71 



1 + ^a 9 (s) ) C eS 

71 



10 • 



(2.5) 



(2.6) 



(2.7) 
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The Cf ew have the advantage of encoding all dominant matrix-element corrections, which lead to 
an explicit s dependence in all of them. 

The functions Oi and 7i(s) were calculated in [11120]; their analytical formulae are given in an 
appendix. The terms (7j(s) take into account universal O (a s ) bremsstrahlung, and the correspond- 
ing infrared (IR) virtual corrections proportional to the tree- level matrix elements of C^-io- The 
remaining (finite) non-universal bremsstrahlung connected with these operators are embedded in 
rate and FB asymmetry through Tj(s). The additional finite terms Tj(s) are rather small, especially 
for large values of s (|xj(s) | < 0.5 for s > 0.3). The finite bremsstrahlung corrections, not related 
to 07_io ® 07-io, are encoded in £r,fb(s) and are substantially smaller. A complete evaluation 
of Sr(s) can be found in [TH] . where its effect is shown to be at the 0(1%) level. The effect of 
Sfb(s) is shown to be below 1% [21]. The coefficients Cj S : w , including the one-loop matrix-element 
contributions of 0i_6, are defined in close analogy with those in Ref. ^1] and are written down in 
an appendix as a function of the true Wilson coefficients Oj. Finally, explicit expressions for the 
latter, evolved down to the low-energy scale, can be found in [T7j . 

The other explicit 0(a s ) terms in (|2.7j) are due to virtual corrections that are infrared-safe. In 
particular, the two-loop functions F^'^ and the one-loop functions Fjp" 9 ' are the matrix elements 
of 0\ ; 2 and Og, respectively, including first order a s virtual corrections. These functions have been 
computed in Ref. [B] for small s. As for arbitrary dilepton mass, these functions will be a part of 
the main fare in this article and will be discussed at great length for the first time in subsequent 
sections. 

2.4 Organization of the perturbative ordering 

As mentioned earlier, the so-called LL order as conventionally labeled is not well justified numeri- 
cally, since the formally leading 0(1/ a s ) term in Cg ew is much closer to being an 0(1) term. For this 
reason, it has been proposed in Ref. to use a different counting rule, where the 0(l/a s ) term of 
Og ew is treated as 0(1). We also subscribe to this approach. Although not consistently extendable 
to higher orders, it is operationally well defined for the present status of the calculation. Within this 
approach, the three Cf ew and the two observables [R(s) and Afb(s)] are all homogeneous quantities, 
starting with an 0(1) term. Then all (7j, and 5i functions are required for a NNLL analysis of 
both R(s) and Afb(s). On the other hand, the missing two- loop matrix element of the operator 
Og and the three- loop mixing of the four-quark operators into Og ^H], which otherwise should be 
a NNLL contributor to the dilepton mass spectrum, are now formally of higher order within this 
new counting. 

3 Calculation of the two- loop matrix elements of Q\,Gi 
3.1 Two- loop diagrams 

The relevant two-loop Feynman diagrams with 0\ and O2 insertions are shown in Fig. ^ They 
will be organized into five gauge-invariant subsets. This is useful because gauge cancellations occur 
within each subgroup, and gauge invariance for each subset is a useful check of the calculation. 
We note that there is another subgroup of two-loop diagrams, shown in Fig. |21 However, only the 
diagram where the virtual photon is attached to the charm loop is non-vanishing. The latter two- 
loop diagram factorizes into two one-loop diagrams and can be included in the calculation of the 
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virtual and bremsstrahlung contribution of the operator Og due to their sharing the same Lorentz 
structure. 

This is most effectively done by the redefinition of the Wilson coefficient (see IA.4I in the ap- 
pendix). When calculating the corrections to the operator Og using these modified Wilson coeffi- 
cients, the contributions of the diagrams given in|2]are automatically taken into account. We also 
note that by gauge invariance the QCD corrections to the matrix elements of the operators 0\^ 
can be written as 

(s£+r \O t \b) = -f- (F?\6g) tree + FP(& 7 ) tree ) . (3.1) 



3.2 Method 



Within the B — > X s 7 calculation at NLL, the two-loop matrix elements of the four-quark operators 
0\ and O2 for an on-shell photon were calculated in ^3] using Mellin-Barnes techniques. This 
calculation was extended in 6\ to the case of an off-shell photon with the help of a double Mellin- 
Barnes representation. We are reminded that these matrix elements form an integral part of the 
NNLL analysis of the decay B — > X s t r i~ . The double expansion is in the dilepton mass s = q 2 /rrib2 
and the mass ratio m 2 c /m 2 . Thus, the validity of the analytical results given in 6\ is restricted to 
small dilepton masses s < 0.25, because cc thresholds will be crossed beyond that. 

We follow here a different strategy to calculate these two-loop matrix elements for arbitrary 
dilepton mass in our present NNLL work. 3 We use a universal method, which can evaluate any 
two-loop diagrams of general external kinematics and internal masses semi-numerically. For its 
implementation, the diagrams are processed with a computer algebra program. The aim of the 
various algebraic manipulations to follow is to render the diagrams to a standard form, which will 
be further integrated numerically in the second stage of the analysis. We used two independent 
versions written in FORM and in Schoonship, which provide a powerful check on the algebra and 
consistency. Let us describe the individual steps of the calculation in more detail. 

First, all two-loop diagrams are converted into sums of sun set type integrals and their mass 
derivatives: 

{{p + k) 2 + m\) ai (q 2 + ml) 012 (r 2 + ml)" 3 

where r = p + q, and p and q are the two independent internal momenta. This is done by the 
use of Feynman parameters {X} and appropriate shifts in the variables p and q. The effective 
masses m\ 2 3 and the effective momentum k are polynomial functions of physical masses, external 
kinematics, and Feynman parameters associated with the diagrams (see Fig. |3J). By using mass 
derivatives, there is in principle only one basic set arising from [oti = a 2 = = 1] we need to know, 
although we shall instead use [a± = 2, «2 = «3 = 1] as the basic set. The reason for the latter choice 
is that this set is less singular and it allows for neater integral representation of the corresponding 
scalar integrals, which makes them more suitable for numerical evaluation. 

The second analytic step includes Lorentz decomposition of the tensor structures and isolation 
of the scalar integrals, and use of differential recursion relations to reduce the scalar functions to a 
set of ten master scalar functions. The tensor reduction is done by decomposing the loop momenta 



d>d n 1 ^ ■ Lw? , Jwo , -lw (3-2) 



3 As mentioned in the introduction, the high-s region above the cc resonances is experimentally also an important 
kinematic window since the efficiency is high there; thus, a comparable number of events will be collected there as 
in the low-s region. However, as we will discuss in sectional one encounters larger non-perturbative corrections in 
this region. 
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a) 



b) 




e) 

Figure 1: Two- loop Feynman diagrams relevant for the virtual QCD corrections corresponding to the operators 0\ 
and 02- They can be organized in five gauge invariant subsets. 
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Figure 2: Two-loop Feynman diagrams relevant to the virtual QCD corrections corresponding to the operators 0\ 
and 02, which can effectively be taken into account by a redefinition of the Wilson coefficients. 




Figure 3: Expressing generic massive two-loop Feynman diagrams as integrals over sunset-type functions. 



10 



p and q in the numerator into components parallel and orthogonal to the external momentum k, 



A 1 



p ■ k 



9 



(3.3) 



Tensor integrals with an odd number of transverse loop momenta p± and q± vanish, while the even 
ones will be contracted basically with the metric tensor for further simplification. The resulting 
scalar coefficients of the tensor decomposition are integrals of the following form: 



Pi^mlmlml-k 2 ) 



d n pd n q 



(p-k) a (q-k) 1 



[(p + k)2 + m 2 ] 2 (q 2 + m?,)((p + q) 2 + ml 



(3.4) 



where a + b < 3 in renormalizable theories; tensor integrals with more than three Lorentz indexes 
can be derived in a similar fashion but will not be needed. In general, one can prove [2E] that any 
two-loop diagram in renormalizable theories can be decomposed by this algorithm into an expression 
involving only ten scalar integrals, denoted by Ti.i, i — 1 • • • 10 in the following. They are linear 
combinations of integrals of the form ()3.4j) with a + b < 3. Let us reiterate that given enough 
computing power and up to possible infrared issues, which we shall mention below, the algorithm 
is applicable to any two-loop process. 

After this algebra step, the ten scalar integrals are integrated over p and q. An important feature 
is that the UV poles are isolated and one arrives at one- dimensional finite integral representations 
over a variable < x < 1 of four elementary functions of rn 2 23 -,k 2 and x, which makes their 
numerical evaluation highly efficient and precise. As anticipated above, for this latter step our 
special choice ol\ = 2, ai = = 1 as a basis of integrals is crucial; with this choice the ten basic 
scalar integrals TCi are logarithmically divergent in the ultraviolet and these UV divergences are 
distinctly separate and manifestly exposed. The remaining finite parts, denoted by hi, possess one- 
dimensional integral representations and will be displayed explicitly in an appendix. The kernels of 
hi are, interestingly enough, moments in the variable x of four elementary functions: 



g(m 1 ,m 2 ,m 3 ;k 2 ;x) 



fi(m 1 ,m 2 ,m 3 ; k 2 ; x) 



f 2 (mi,m 2 ,m 3 ; k 2 ;x) 



/ 3 (mi,m 2 ,m 3 ; k 2 ; x) 
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Once the analytical procedure is done as described, each original Feynman diagram is expressed 
as an integral over the set of Feynman parameters {X} introduced earlier. The integrand itself 
consists of a sum of the special functions hi (which are themselves one-dimensional integrals of 
the four elementary functions ()3.5Jl ) and possibly also of some trivial functions such as logarithms 
and rational functions of the kinematical invariants Tnf 2 3> and k 2 . Within our method, all these 
integrations generally are to be performed numerically. We shall not dwell on the details here. 
Suffice it to say that the analytic structure of the functions hi is well understood, so that the 
integration paths {C} can be moved into complex planes to effect better numerical convergence 
and, more importantly, to yield amplitudes on the physical sheet. Because we are interested in a 
high accuracy and efficiency routine, we have used an adaptive deterministic integration algorithm. 
Such integration routines are very accurate, provided that the integrand is smooth enough and that 
the dimensionality of the integral is not too large. The integrand itself is, of course, an analytic 
function along any properly chosen complex integration path of {C}, and therefore in order to 
preserve this smoothness, it is advantageous to optimize the choice for a smooth integration path in 
the numerical work as well. One should be made aware that the integrals over Feynman parameters 
must be performed along a complex integration path that is consistent with the causality condition. 
This path is computed automatically by using spline functions such that both the path itself and 
its Jacobian are smooth functions. Moreover, we note that, in the problem at hand, we shall be 
dealing with three- fold numerical integrations at most. 

We would like to mention that we perform minimal subtractions on all divergent subgraphs. We 
have checked that the anomalous dimensions so obtained for the operator mixing matrix elements 
agree with what is known in the literature. This will be further elaborated in a subsequent section. 
We have explicitly checked that each group of diagrams is gauge-invariant and have been further 
reassured by their numerical stability. 

In the case of the calculation of the two-loop matrix elements of B — > X s £ + £~ decay, we deal with 
three kinematic variables: the charm mass, the dilepton invariant mass, and the subtraction scale. 
In order for the result to be usable for phenomenological studies, in particular to be implementable 
into a Monte Carlo simulation, we need to cover this three-dimensional kinematic space. A real-time 
calculation of the two-loop matrix element by numerical integration is far too slow to be implemented 
directly into a Monte Carlo simulation of the experimental set-up. With the present-day processors, 
an alternative of calculating in advance a comprehensive grid of integration points that cover the 
whole three-dimensional kinematic range, storing them, and then interpolating between them, seems 
to be the most efficient way because of the semi-numerical nature of our whole approach. All these 
considerations led us into writing a program that calculates the two-loop matrix elements efficiently 
and accurately, and which is fast enough to be incorporated into a Monte Carlo simulation. We 
selected a grid of 38 x 3 x 3 integration points for both the electric and the magnetic components 
of the two-loop virtual corrections. Each of the 684 integration point yielding values for the form 
factors was calculated with a relative precision of 10~ 3 . We used the CERN Linux cluster to perform 
this calculation, and the CPU usage was approximately 3 days on 33 processors (mostly 850 MHz) 
running in parallel. 

Finally, we would like to mention a caveat in the semi-numerical algorithm presented here. In 
a general process, infrared singularities will most likely occur and will lead to infrared divergences 
in our integral representations. In such cases, it is more efficient to first separate out the infrared 
parts of the two-loop diagrams in an analytically manageable form. Fortunately, for the process 
at hand all relevant diagrams (a)-(e) in Fig. are infrared-finite; therefore the algorithm is most 
suitable for this specific application. 
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3.3 Unrenormalized results 



In the following we show plots of our results for the finite (e°) parts of the unrenormalized (naked) 
two-loop Feynman diagrams of Fig. which give matrix elements to the operators 0\ and O2 within 
the MS scheme. The finite counterterm contributions are not included here but will be discussed 
in the next subsection. Our main purpose is to compare our results with those of Ref. jB] where 
Mellin-Barnes expansion techniques were used. 

Each plot in Fig. 0] represents one of the five gauge-invariant diagram subsets given in Fig. ^ 
We want to bring attention to the complete range of the dilepton mass spectrum s = [0, 1] in all 
the plots, which is a salient point of this discussion. For each we plot the electric (F^) and the 
magnetic (F^) contributions, defined in ()3.1)) separately in the left and right columns, except for 
the subset (e), which has no contribution in the first column. The calculation of Asatrian et al. is 
shown as successive approximations, showing that it converges toward our result. Their results are 
valid only under the cc threshold, which we demarcate by vertical lines in the figures. Generally, 
the real part of our results is given by a solid line and the imaginary part by a dashed line. All the 
other lines are successive approximations in the momentum expansion series of Asatrian et al.. 

In Fig. El we zoom in to magnify the way the momentum expansion converges toward our exact 
numerical solution within the low-s region. Here the stars are the actual points obtained from our 
numerical output, which are joined by the continuous lines obtained from the interpolation program. 

The comparison can be summarized as follows: there is good agreement between our results for 
each diagram set and the double expansion of Asatrian et al. in the region below the cc threshold. 
This agreement within the low-s region provides a strong confirmation of our numerical method. One 
notices that, as a general rule, the less singular the threshold behaviour of the diagram is, the better 
the momentum expansion converges toward our exact numerical result. For subsets (a) and (b) the 
expansion converges best because of the lack of a nontrivial threshold. For (c) and (d) the threshold 
is relatively mild, and the convergence is intermediary. For the gauge-invariant subset (e), the 
threshold is quite singular and thus we notice the poorest convergence of the momentum expansion. 
This singular behaviour is mostly due to the charm self-energy-type diagrams, as evidenced by a 
much milder disagreement of the two methods after the mass counterterm is added. This also makes 
the agreement between our final physical result and the momentum expansion result better than 
what can be inferred from Fig. |3] alone. The actual agreement of the two calculations within the 
low-s region is compatible with the error of our numerical integration accuracy. We may find it 
interesting that, for the gauge-invariant subsets (a) and (b) the expanded results of [U] are actually 
valid beyond the cc threshold, perhaps because these diagrams have no threshold cuts in the way. 

Before leaving this subsection, we want to remark that by our method we have reproduced 
numerically the values given by the analytical results on the two-loop matrix elements of Oip in 
the B — > X s 7 mode presented in [T3], with an accuracy well below 1%. 

3.4 Counterterm contributions 

The counterterms to the matrix elements of the operators 0\ and O2 will serve two purposes. They 
give renormalization to the QCD parameters and they account for operator renormalization due to 
mixing. Although the counterterms expanded in the variable s can be found in [U], for us we need to 
generalize to arbitrary dilepton mass. Before giving our new final analytical results, we commence 
with a short discussion of the various bits and pieces that must be put together. We follow the 
conventions in jB], including the convention d = 4 — 2 e and their convention regarding the evenescent 
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Figure 4: Plots of the UV- finite part of the Feynman diagram subsets shown in Fig.^ The two columns correspond 
to the electric and the magnetic form factors; see (|3 .11) . There is no electric contribution from the diagrams 1 e). We 
plot the real (solid line) and the imaginary (dashed line) parts of our exact numerical integration result, along with 
successive approximations in the momentum expansion series of the result of Asatrian et al. which is shown only 
below the cc threshold. For more details see main text. 
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Figure 5: Convergence pattern of the momentum expansion solution of Asatrian et al. toward our numerical 
integration solution. The stars denote actual integration points, and the solid line connecting them is an interpolation. 
For more details see main text. 
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operators. The anomalous dimensions describing the operator mixing are known (see ^7J ITHj ) 
and are given in the appendix. We take this opportunity to point out a very useful feature of our 
semi-numerical method in that the UV divergences are separated and given analytically through 
the whole calculation. Because of this, we were able to check the two-loop mixing results explicitly. 

The list of non-trivial counterterm contributions to the functions and in Ijfl.lj) . which 
are add to those from the naked diagrams given in Fig. ^ is the following: 

• The two operators 0\ and O2, which non-trivially mix into the four-quark operators at the 

(7) 

one-loop level, inducing additional counterterm contributions. We denote them by F^_J Aquark 
and F^ Aquark . They are given by 

H ^~a>ij( S t + £~\Oj\b) one-loop = -(^( F i^4quark(^7)tree + F-^ 4quark {0 9 ) tree)) , (3.6) 
3 

in which j runs over all four-quark operators. This expression instructs us to calculate the 
one-loop matrix elements of the four-quark operators up to the e 1 terms for arbitrary dilepton 
mass. 

• The analogous one- loop mixing of the operators 0\ and O-i into Og. They are from the 
two-loop diagrams shown in Fig. |21 



The one and two- loop mixing into Og, which is connected with the renormalization of the 
explicit o 
leads to < 
given by 



explicit coupling constant in the definition of Og. They appear in the diagrams of Fig.Q] This 
leads to an additional contribution to F^ (but not F- 7 ') in (|3.1|) . which we call F^ 9 . It is 



e e 

where we have used the coupling renormalization 



„22 ,,12 _11 a 
ri (9) ( a i9 , a i9 \ a i9 PO / o *y\ 

fii^g = -(— + —) (3.7) 



Z gs — 1 — — ^ ^ — h 0{a 2 s ) with (3q = 11 — — for five active flavours. (3.8) 
An 2 e 3 

• Then, there are QCD mass counterterm contributions from the renormalization of the charm 
mass within the matrix elements of 0\ and O2, which we denote by F^ c and . They 
are most easily given by replacing the charm mass m c by Z mc x m c in the one-loop matrix 
elements of 0\ and 02- We use the pole mass renormalization in our calculation: 

z ^-t:I(W^ +4 ) + °^ (3J) 

In the following, we give the complete finite (e°) counterterm contributions from the two- loop 
diagrams in Fig.^ We write down our new final results for the relevant finite (e°) parts for arbitrary 
dilepton mass: 

F$ = F(% + F}% ark + F^ c , k = 7,9,2 = 1,2. (3.10) 
F $t U° = ^ + ^B(s) (3.11) 



51 81 
27 M 27" 



Fgeileo = -^-^B( S ) (3.12) 
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(3.14) 



The functions 5(a), B x (a), B xx (a), -82(a), B2 X (cl), B2 XX (a) and also the functions C(a), (^(a), 
^2(0), C*2x(a), C2xx(a) are defined in the appendix. We have also introduced the variable 
s c = q 2 jm\ and we have defined Z M = ln(m 2 //z 2 ), Z c = ln(m 2 /m 2 ), and Z s = ln(s). 

We have explicitly checked that all the e~ 2 and e _1 terms coincide with the results given in jB]. 



4 Calculation of the 0(a s ) virtual corrections to the 
matrix element of Og 

We present here a short description of our new calculation of the matrix element of the operator 
Og for general s. Besides the contributions from the naked diagrams shown in Fig. El there is also 
a counterterm contribution due to the mixing of Og into 7 : 

n> 16 

(s£ + £-\SZ 87 8 \b), 5Z 87 = --^-. (4.15) 

47r 9e 

The functions F 8 7 ^ and F 8 9 ^ in ()2.7j) are then defined by the renormalized matrix element of 8 : 

(se + £-\C 8 8 \b) = c£\-^)(F 8 {9 \d 9 ) tree + FP(6 7 ) tree ), (4.16) 
Keeping the full s dependence, we find 
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Figure 6: One-loop Feynman diagrams relevant to the virtual QCD corrections corresponding to the operators Os- 
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(4.18) 



If we expand our results for small s, we recover the results given in j^]. 



5 Non-perturbative contributions 

5.1 Generalities 

Non-perturbative contributions in B — > X s £ + £~ transitions can be divided into two main categories: 

• Aqco/^b corrections in the relation between the partonic b — ► s£ + £~ amplitude and inclusive 
hadronic distributions; 

• non-perturbative effects associated with the cc intermediate state: B — ► X s cc — > 

The heavy-quark expansion, which led us to evaluate the first type of contributions, is rapidly con- 
vergent and leads to small corrections for sufficiently inclusive observables. A consistent treatment 
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of the second type of effects requires to impose kinematical cuts to avoid the large non-perturbative 
background of the narrow cc resonances. These two requirements are somehow in conflict. As a 
result, we can perform reliable predictions of B — > X s £ + £~ transitions, both in the low- and in the 
high-g 2 regions, but magnitude and error of the non-perturbative corrections are enhanced with 
respect to their natural size. 

The enhancement of corrections is particularly sizeable in the high-g 2 region, because 

of two main drawbacks: 

• the l/m b expansion breaks down in the limit g 2 — > ml [2UI2Z]; 

• the q 2 cut introduces a sizeable effective correction linear in l/ra&, through the relation between 
hadronic and partonic phase spaces. 

The first problem implies that in the high-g 2 region the differential distribution in q 2 cannot be 
predicted in perturbation theory. This non-perturbative distribution has nothing to do with the 
so-called shape function, or the kinetic energy distribution of the heavy quark inside the hadron j2Zj . 
However, similarly to the latter, the q 2 distribution near the end point is a non-perturbative function, 
which must be determined from data. What can still be predicted with reasonable accuracy, for a 
sufficiently low cut-off g^ in , is the q 2 > g^ in integral (or the full inclusive distribution for q 2 > g^ in ). 

The second drawback is common to all observables that require kinematical cuts (in practice to 
any experimentally accessible inclusive observable in B decays). Since any kinematical cut on the 
final state must be expressed in terms of the hadron mass Mb, we cannot avoid the linear l/m& 
corrections that arises from the relation 



M B = m b 



' , A Ai + 3A 2 / Aqcd 

1 + 2 + O — , a s 

nib ^ m b \ m b 

2 



(5.19) 



This problem is substantially enhanced in the high-g region because of the smallness of the available 



phase space: here the relative correction between the hadronic phase space [~ (Mb — v<Zmiii)] an d 



the partonic one [~ (m b — yc^m)] becomes an 0(1) effect. 

As we shall discuss in detail in the following, these two drawbacks fit within a common picture: 
the heavy-mass expansion in the high-g 2 region is an effective expansion in inverse powers of 



™heavy = m b X (1 - y/s~) , (5.20) 

rather than m b . This expansion is justified, but it converges less rapidly than the usual series in 
Aq C d/w 6 . 



5.2 Aq CD /?77| and Aq CD /?77| corrections 

The Aqcd/^6 corrections can be systematically investigated in the framework of the heavy-quark 
expansion and, in particular, by means of the heavy-quark effective theory (HQET) [2E]- The two 
main distributions, R(s) and A FB (s), are not affected by linear corrections and the leading effects 
of 0(AQ CD /m 2 ) can be described in terms of the two expectation values 

_ {B\h(tD) 2 h\B} _ 1 (B\hga ■ Gh\B) _ M|, — M| 

Al " 2Mb ' 2 ~ 6 2Mb ~ 4 ' (5 } 
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where h is the heavy-quark field in the effective theory. The explicit expression of these corrections, 
which have been computed in Refs. [2HI2Z1 (see also Ref. (23), is 
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(5.22) 



+2Re [C%r(s)*ar(s)] (7 + 10s - 9s 2 ) + ^^A FB (s) 



+ 3mJJl^y AMs) 



(5.23) 



In both cases we have taken into account also the 1/m 2 terms arising from the semileptonic nor- 
malization, namely 
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(5.24) 
(5.25) 



This normalization is responsible for the absence of explicit 0(Aqcd/^&) corrections, since it cancels 
any explicit dependence from m&, and it is also responsible for the absence of any dependence from 
the kinetic energy of the 6-quark (~ Ai) in R(s). 

As can be seen in Fig. the relative corrections 5 1 / m 2R(s)/R(s) and S 1 / m 2AFB(s)/ApB(s) 

* b • b 

are of the order of a few percent in the small s-region, apart from the obvious divergence in 
5 1 / m 2A FB (s) / A FB (s) due to the zero of A FB (s) (the shift in the position of the zero amounts only to 
an increase of about 0.9%). However, in both cases the non-perturbative corrections become quite 
sizeable in the large s-region and the 1/mb expansion breaks down close to the s — > 1 endpoint. 
The nature of this singularity has been discussed in detail in Ref. [2Z] • As usual, the HQET cannot 
be applied in corners of the phase space of 0(AQ CD /m{>), where the kinematics forces the final 
hadronic state to assume soft configurations. However, this particular case is rather different from 
the well-known examples of the photon-energy endpoint in B — > X s r y 1 or the lepton-energy endpoint 
in B — > X c iv. There only the hadronic invariant mass is constrained to be semi-soft (k 2 ~ mtAQco) 
and the breakdown of the HQET is cured by means of a resummation of singular terms which leads 
to the shape function, or the universal non-perturbative distribution of the 6-quark kinetic energy 
inside the hadron. Contrary to these examples, the kinematical constraint corresponding to the 
dilepton invariant-mass end-point, namely m b — ^fq 1 = O(Aqcd), forces the hadronic system to 
have both soft momentum (k ~ Aqcd) and soft (k 2 ~ Aq CD ) invariant mass. This implies that 
no resummation can be applied and that the singularity has nothing to do with the kinetic-energy 
distribution of the 6-quark [27j . 

The fact that the 0(Aqcd/^&) corrections are not under control for s — > 1, does not prevent us 
from performing reliable predictions for the partially integrated branching ratio (and FB asymmetry) 
in the high-s region, provided we choose a sufficiently low cut-off s m ; n . The main issue is which is 
the maximal allowed value for Smin. 
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Figure 7: Relative corrections due 0(Aq CD /m^) effects: 5 1 / m 2R(s)/ R(s) (dotted) and Sy m a A F b(s)/A F b(s) 
(dashed). 



A 


Af (GeV 2 ) 


Af (GeV 2 ) 


Pi (GeV 3 ) 


0.40 ±0.10 


-0.15 ±0.10 


0.12 ±0.02 


0.06 ±0.06 



Table 1: Input values of the HQET parameters used to estimate linear, quadratic and cubic corrections in the 1 /nib 
expansion. 



Once we impose a constraint of the type s > s m j n , the inclusive sum on the hadronic final state 
is limited to systems with invariant mass up to the effective scale m^ vy in (I5.20|) . For this reason, 
in this partially integrated observables we expect an effective expansion ruled by inverse powers of 
m heavy> rather than nib. This naive expectation is confirmed by the detailed analysis of Ref. "3*0*] . 
applied B — > X u lv decays. There the dilepton-invariant-mass cut necessary to avoid the B — > X c iv 
background leads to an effective expansion in inverse powers of m c , rather than nib j30j . 

Because of these general arguments, we expect that the expansion should still be reliable, al- 
though with a slower convergence, for s min > 0.6 (corresponding to ^f avy ~ 1 GeV). To address 
this issue in a more quantitative way, we shall look at the explicit expression of the Aqcd/ 777 ! cor- 
rections in R(s) jSI]. At this order we need to introduce seven new hadronic matrix elements. Five 
of them lead only to a redefinition of the AQ CD /m 2 couplings Ai^- In particular, the contributions 
proportional to T\..a and p 2 , in the notation of Ref. $2\, are obtained from Eqs. (|5.22j) and (|5.23j) . 
with the replacement pH] 13*2*] : 

Ax - Af = Ax + HL±^, Aa-Af = A 2 ±^±^-^. (5.26) 
nib mb mb 

These terms do not spoil the convergence of the I /nib expansions, independently of the s m i n cut, 
provided the naive chiral- counting expectation 7i... 4 ~ pi 2 ~ Aq CD is respected. More delicate is 
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the issue of the contributions proportional to p\ and j\ [HI] : 
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-3s 3 + 17s 2 - s - 



+ 8A(/ 1 )5(l-s) 



4ia n 



+ 8A(A)(J(l-s) 



|C 9 new ( 



+ ic 
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where 



2(1 -s) 

77 - 88^ + 24z 2 



+ 8A(A)5(l-s) 



Re [C 7 new (s)C 9 new (s) 



(5.27) 



^(s) = 77 - 88^ + 24^ - 8z s + 5^ 4 + 48 In z + 36z 2 In z (5.28) 

arises from the semileptonic normalization and A(/i) is a local contribution that cures the singu- 
larity of jQds'R(s') for s — > l. 4 When integrated in the high-s region, Eq. ()5.27j) leads to a huge 
coefficient for the pi/mf correction, much larger than the already sizeable A2/m 2 term [HI]. How- 
ever, everything looks very reasonable, once we introduce the effective scale in (|5.20j) . For s m ; n ~ 0.6 
we find 



ds R(s) 



1.6A 9 



m 2 (l - y/T r _ 



+ 



1.8 Pl + L7A 

m f(l ^/^min, 



X 



f 1 ds R(s) 



Img,— >oo ' 



(5.29) 



which perfectly confirms our expectation of an effective expansion in inverse powers of w|f avy - Ac- 
cording to the input values in Tabled which are consistent with recent experimental determinations 
(see e.g. Ref. jHH]), and setting f x = 0, 5 the numerical size of the term between square brackets in 
flo"^]) is [1 -0.08 ±0.08]. 



5.3 The A/rrib correction 

As anticipated, even though R(s) and Afb(s) are not explicitly affected by linear corrections in the 
l/m& expansion, the physical observables defined in terms of a q 2 cut are sensitive to the A/rrib term 
via the relation (|5.19|) . This term - or equivalently the uncertainty on the value of m;, - represent 
at present the largest source of non-perturbative uncertainty in the high-g 2 region. Choosing as 

= 0.6, the physical observable defined in terms of q 2 in can be written 



reference cut the value s r 



as 



Rcutil; 



2 - ) 

mm/ 



dq 



2 dT(B 



T(B -> X r ev) 



6.2 



^min 

ml 



-0.6 



O 



mi 



-0.6 



0.6 



ds R{s) 



which implies 



SR 



cut 



Rr 



7.4 



Srrib 
m b 



(5.30) 



(5.31) 



4 The cut-off-dependent coupling A(/i) is related to the cut-off-independent parameters /1 and pi, defined as in 
Ref [21], by the relation J* ds[l/(l - s) + A(/i)<5(l - s)] = -fi/pi. 

5 To fix this quantity requires more restricting information, but we assume its contribution is within the present 
uncertainty. 
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This means that an error Srrib = 0.1 GeV (corresponding to the uncertainty on A in Table HJ, leads 
to a ~ 15% error on R cn t- 

5.4 I /ml non-pert urbative corrections 

The second class of non-perturbative effects relevant in B —>■ X s £ + £~ decays are the long-distance 
corrections related to cc intermediate states. These originate from the non-perturbative interactions 
of the cc pair in the process B —>■ X s cc X s £ + £~ . If the dilepton invariant mass is near the first 
two J PC = 1 cc resonances and this effect is very large and shows up as a peak in R(s). 
However, one can easily eliminate this background by suitable kinematical cuts. More delicate is 
the estimate of the long-distance effects away from the resonance peaks. 

An interesting approach which avoid double- counting problems is the one proposed in Ref. |33| 
(KS approach). Here, in order to take into account charm rescattering, the correction to Cg induced 
by b —>■ ccs operators is estimated by means of experimental data on a(e + e~ — > cc hadrons) using 
a dispersion relation. To be more specific, the function h(z, s) appearing in ()A.5j) is replaced by 

h(z, s) — h(z, 0) + -P £ ds'^-j^ + ^-R~ d (s) , (5.32) 

where R^ ad (s) = a(e + e~ — > cc)/cr(e + e~ — > /i + yU~) and s c is the cc threshold. This method is 
exact only in the limit where the B — > X s cc transition can be factorized into the product of sb 
and cc colour-singlet currents (i.e. non-factorizable effects are not included). The non-perturbative 
corrections estimated using this approach are extremely small in the perturbative windows s < 0.25 
and s > 0.6 [SHI- Fo r the integrated branching ratios one finds an increase of 1 — 2% in the low-s 
region, while the effect in the high-s region is far below the uncertainty of the corrections. 

A systematic and model-independent way to estimate the non-factorizable cc long-distance ef- 
fects far from the resonance region is obtained by means of an expansion in inverse powers of the 
charm-quark mass (3HIEZI- This approach, originally proposed in j^Hl to evaluate similar effects 
in B — > Xsj decays, has the advantage of dealing only with partonic degrees of freedom. In this 
framework the leading non-perturbative corrections to R(s) and Afb turn out to be 0(AQ CD /m^); 
their explicit expressions can be found in the appendix. Since the factorizable corrections vanish for 
s — >• 0, the 0(Aq CD /ml) effect is expected to be the dominant long-distance contribution for small 
values of the dilepton invariant mass. In this region the relative magnitude is very small (at the 1 
or 2% level) and opposite in sign to the factorizable KS correction. The 0(AQ CD /m^) calculation 
should be reliable also above the resonance region (s > 0.6), where the effect is again very small. 
Similar comments apply to the long-distance corrections for Afb(s). 

6 Phenomenological analysis 

6.1 Branching ratio and dilepton invariant-mass spectrum 

The final results of this work concerning the dilepton invariant-mass spectrum is summarized in 
Fig. |H1 In the upper plot we compare our un-expanded result, without any non-perturbative cor- 
rection, to the expanded result of Ref. U . As can be noted, the expanded result provides a perfect 
approximation to the full calculation up to about the cc threshold. This is of course a good cross- 
check of Ref. [Hj and an important test of our method, which turns out to be essential to provide 
a reliable prediction in the high-g 2 region. The Fig. |H1 also shows that the expanded result differs 



23 



m b = (4.9 ±0.1) GeV m c /m b = 0.29 ± 0.02 



5.0l|j) GeV 



a s {M z ) = 0.119 


«em = 1/128 


\V t * b V ts /V cb \ = 0.97 


m t (m t ) = 167 GeV 



Table 2: Main input values used in the numerical analysis. 



significantly from the exact NNLL result above the threshold as expected. Note that the scale de- 
pendence in the high-g 2 region is very small, therefore our NNLL result provides an excellent level 
of accuracy for the pure partonic calculation. The lower plot in Fig. |H] provides an illustration of 
the non-perturbative effects induced by cc intermediate states (evaluated using the KS approach, 
see Sect. 15. 4|) . As can be noted, the point-by-point corrections in q 2 are quite sizeable in the high-g 2 
window. As discussed in the previous section, only the g 2 -integral can be predicted reliably in this 
region. 

Before discussing the numerical predictions for the integrated branching ratios, we wish to 
emphasize that low- and high-g 2 regions have complementary virtues and disadvantages. Taking 
into account the discussion in the previous section, we can summarize the main points as follows: 

• Virtues of the low-q 2 region: reliable q 2 spectrum; small l/m b corrections; sensitivity to the 
interference of C-j and Cg] high rate. 

• Disadvantages of the low-q 2 region: difficult to perform a fully inclusive measurement (severe 
cuts on the dilepton energy and/or the hadronic invariant mass); long-distance effects due to 
processes of the type B —>■ ^X s — ■> X s + X'£ + £~ not fully under control; non-negligible scale 
and m c dependence. 

• Virtues of the high-q 2 region: negligible scale and m c dependence due to the strong sensitivity 
to |Cio| 2 ; easier to perform a fully inclusive measurement (small hadronic invariant mass); 
negligible long-distance effects of the type B — > ^/X s — > X s + X'£ + £~ . 

• Disadvantages of the high-q 2 region: q 2 spectrum not reliable; sizeable l/m b corrections; low 
rate. 

Given this situation, we believe that future experiments should try to measure the branching ratios 
in both regions and report separately the two results. These two measurements are indeed affected 
by different systematic uncertainties (of theoretical nature) and provide different short- distance 
information. 

In order to obtain theoretical predictions which can be confronted with experiments, it is nec- 
essary to define the two regions with appropriate cuts in q 2 (and not in the partonic variable s, as 
done in most of the previous literature). Concerning the low-g 2 window, we propose as reference 
interval the range q 2 G [1, 6] GeV 2 . The lower bound on q 2 is not essential, but it proposed in order 
to cut a region where there is no new information with respect to B — > X^ and where we cannot 
trivially combine electron and muon modes. The higher cut is essential to decrease the uncertainty 
associated to the cc threshold. 

Taking into account the input values in Table |2l the NNLL prediction within the SM for this 
low-g 2 window is: 

cut Ji GeV 2 T{B -> X c ev) 
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Figure 8: NNLL predictions of dB(B — > X s l + l~)/dq 2 . Upper plot: pure partonic result with full m c dependence 
computed in this work for n=5 GeV (full line), /x=2.5 and 10 GeV (dashed lines); g 2 /(4m^)-expanded result by 
Asatrian et al. UJ, for fi=5 GeV, extrapolated to the full q 2 range (dotted line); the dashed vertical lines indicate the 
positions of the first two narrow ^ resonances. Lower plot: partonic result with full m c dependence for /i=5 GeV 
with (dotted line) or without (full line) factorizable cc corrections computed in the KS approach (see Sect. I5.4|> . All 
other inputs are fixed to the central values in Tabled 
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1 ± 8%| r- ± 6.5%|^ ± 2%| mc ± 3%| mb(cuts) + (4.5 ± 2)%| 1/m g - (1.5 ± 3)%| ( 

= (1.52 ± 0.18) x 10" 5 . (6.33) 

Between square brackets we have reported all the uncertainties and non-perturbative corrections 
discussed in this work, evaluated according to the input values in Table d and |2] The error denoted 
by r s i corresponds to the theoretical uncertainty implied by the T(B — > X c ev) normalization which, 
in turn, is dominated by the uncertainty on m c . In principle, alternative normalizations such as the 
one proposed in Ref. jHU] could be used to reduce this error. In any case, this uncertainty should 
be regarded as a parametric error which can be improved by additional independent measurements. 
The small error denoted by m c correspond to the m c -dependence of F(B — > X s £ + £~) only (ignoring 
the normalization): as can be noted, this is almost negligible. 

As already pointed out, our calculation of the matrix elements provides a completely independent 
check of the results of Ref. jH] for the low-g 2 region. In particular, we confirm the reduction of the 
scale uncertainty from ±13% to ±6.5%, in the branching ratio, once these NNLL corrections are 
included (see Fig. EJ). The difference in the central value of R l °^, compared to the result of Ref. 0, 
is entirely due a different definition of this observable and to differences in the input values. 

Concerning the high-g 2 window, we propose as reference cut q 2 > 14.4GeV 2 , which leads to the 
following first NNLL prediction of the high dilepton mass spectrum: 

hi gh _ f , 2 dT(B — > X s ( + 



iff = / -^ = 4.09x10- 

Jq 2 >lAA GeV 2 Y{B — > X c eu) 



(my, — 4.9 GeV\ 

1 ± 8%| rsj ± 3%|„ + 0.15 ( b Q1GeV ) - (8 ± 8)%| 1/m ( 2! a) ± 3% 
(3.76 ± 0.72) x 10" 6 . (6.34) 



Here the explicitly indicated m& dependence induces the largest uncertainty. At present this is 
about 15%. However, significant improvements can be expected in the near future in view of more 
precise data on other inclusive semileptonic distributions. Note that, as anticipated, in this region 
the pure perturbative uncertainties due to scale and m c dependence are very small (the latter has 
not been explicitly indicated being below 1%). The impact of the NNLL corrections computed in 
this work for the high q 2 region is a 13% reduction of the central value and a significant reduction 
of the perturbative scale dependence (from ±13% to ±3%, see Fig. [0J). 

As also shown in Fig. El this reduction brings the central value of the full NNLL prediction very 
close to the partial NNLL result obtained without the F^s, functions, provided in the latter case 
the renormalization scale is set around 2.5 GeV. This observation, which has already been made in 
Ref. j3U] for the low-g 2 window, applies remarkably well also in the high q 2 region. Obviously, our 
full NNLL calculation provides a fundamental ingredient to justify this procedure for the central 
value and, especially, to obtain a clear estimate of the residual scale uncertainty. 

It must be stressed that two non-negligible source of uncertainties have not been explicitly 
included in Eqs. f|6.33|) and ()6.34|) : the error due to m t (and the high-energy QCD matching scale) 
and the error due to a cm (or better the error due to higher-order electroweak and electromagnetic 
effects). The first type of uncertainty has been discussed in detail in [T7j, and it amounts to ss 6%. 
As far as the uncertainty on higher-order electroweak corrections is concerned, the error is also 
expected to be at the level of a few percent, but a consistent estimate of these effects is beyond the 
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Figure 9: Complete NNLL result for dB{B -> X s £+£-)/dq 2 for ^=5 GeV (full line), fi=2.5 and 10 GeV (dashed 
lines), vs. the partial NNLL result obtained by neglecting the F^'g functions computed in this work, for /i = 2.5 GeV 
(lower dotted line), [i = 5 GeV (middle dotted line), and [i = 10 GeV (upper dotted line). 



-> X c ez/) = (10.74 ± 0.24)% @2|, we finally obtain: 

; q 2 E [1,6] GeV 2 ) = (1.63 ± 0.20) x 10" 6 , (6.35) 
; q 2 > 14.4 GeV 2 ) = (4.04 ± 0.78) x 10~ 7 . (6.36) 

At the moment, these two predictions cannot be directly compared with experimental data. To 
facilitate the comparison with the recent results by BELLE and BABAR [7J |H], we present also 
an estimate for the interpolated partonic rate. The latter is defined as the integral of the partonic 
rate (full line in Fig. |HJ) over the full q 2 spectrum (including the resonance region), starting from 

B[B -> X s £ + e~ ; q 2 > 4m 2 ) = (4.6 ± 0.8) x 10~ 6 . (6.37) 

The central value of this estimate, based on our NNLL calculation, takes into account the power 
corrections in the clean perturbative windows, while its error includes a guestimate of the systematic 
uncertainty due to the extrapolation (evaluated using the KS approach). The difference of our 
central value with respect to a similar estimate for the interpolated partonic rate presented in 
Ref. is almost entirely due to parametric differences in the input values (most notably, to the 
choice of a em ) and is well within the theoretical errors given in ([6.37)1 . 

Our estimate in ()6.37|) compares well with the recent experimental world average |41j : 

B[B -> X s £ + O cxp = (6-2 ± l.llli) x 10" 6 . (6.38) 

6 Choosing as reference value for a em the value at the electroweak scale, we should have minimized the impact 
of the dominant electroweak matching corrections. After this work has been completed, a complete analysis of 
higher-order electroweak corrections has been presented in Ref. 



scope of this work. 

Using the world average T(B - 

b(b -> x s e + r 
B(B -> x s i + r 
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Figure 10: NNLL perturbative contributions to the un-normalized FB asymmetry; notations as in the upper plot 
of Fig. E 

In view of data with higher statistics, we stress once more the importance of a future compar- 
isons with the more clean and more interesting predictions in Eqs. (|6.35|) and ([6.35)1 . Contrary to 
Eq. ()6.37|) . the theoretical errors in both (|6.35|) and ()6.35|) could be systematically improved in the 
near future. 



6.2 Forward-backward asymmetry 

The summary plots for the lepton forward-backward asymmetry are shown in Fig. El and ^2 In 
Fig. El we plot the un-normalized differential asymmetry, defined by 



dA 



FBI 



dq 2 



XJ+t 



1 dcQsdi d 2 B(B 
-l dq 2 dcosO 



-sgn(cos 6g 



B(B 



m 



\ m b/ 

while in Fig. we plot the (adimensional) normalized differential asymmetry, defined by 

1 



(6.39) 



a fb (q 



dB(B -> X s £+£-)/dq 2 J-i 



i d 2 B(B - X s 

d cos Ui 



dq 2 dcosOe 



-sgn(cos^) . 



(6.40) 



Most of the comments concerning the errors and the complementary of low- and high-g 2 windows 
discussed before holds also for the forward-backward asymmetry. 

In the low-g 2 region the most interesting observable is not the integral of the asymmetry, which 
is very small due to the change of sign, but the position of the zero. As discussed by several authors 
(see e.g. Ref. [HEDD, this is one of the most precise predictions (and one of the most interesting 
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Figure 11: NNLL perturbative contributions to the normalized FB asymmetry; notations as in the lower plot of 
Fig.! 



SM tests) in rare B decays. Denoting by q 2 the position of the zero, and showing explicitly only 
the uncertainties and non-perturbative effects larger than 0.5%, we find at the NNLL order 



ql 



0.161 x mi x 



l + 0.9%| iy ,±5%| NNNLL 



(3.90 ± 0.25) GeV 2 



(6.41) 



As already pointed out in Ref. jlj, in this case the /i dependence is accidentally small and does 
not provides a conservative estimate of higher-order QCD corrections. The 5% error in 1)6.41)1 has 
been estimated comparing the result within the ordinary LL counting and within the modified 
perturbative ordering proposed in Ref. (see Sect. 12. 4j) . The central value, as well as all the other 
central values in this work, is obtained using the modified ordering of Ref. 

In the high-g 2 window the FB asymmetry does not change sign, therefore its integral represents 
an interesting observable. In order to minimize non-perturbative and normalization uncertainties 
it is more convenient to consider a normalized integrated asymmetry. Applying the same q 2 cut as 
in (J6.34)) . we define 



(^FB)cuts 



q 2 >14A GeV 2 



dq 2 



dB FB (q 2 
dq 2 



q 2 >lAA GeV 2 



dq 



2 dA F B(q 2 ) 



dq 2 



(6.42) 



All parametric and perturbative uncertainties are very small in this observable at the NNLL order 
level. On the other hand, despite a partial cancellation, this ratio is still affected by a considerable 
amount of Aq CD /m 2 and AQ CD /mjj corrections (which represent the by-far dominant source of 
uncertainty). Taking into account the expressions in Sect. 15.21 and separating the contributions of 
the various operators, we find 



(^FB)cuts 



0.42 x [1 - (0.17±0.11) Al 
0.14 ±0.06 . 



(0.42 ± 0.07)a 2 - (0.08 ±0.08; 



(6.43) 
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Appendix 1: Auxiliary functions 

• The function h(z) describing next-to-leading order QCD corrections to the semileptonic decay 
[see Eq. (Q] is given by [23j: 



y _ 1 ~2 , 17 3 

* 3 * " r 3 A 



h(z) = -(1 - z 2 ) (f - ?f z + f z 2 ) + z ln(z) (20 + 90 

+z 2 hi 2 (z) (36 + z 2 ) + (1 - z 2 ) (f - f z + f z 2 ) ln(l - z) 
-4 (1 + 30 + ^ 4 ) ln(z) ln(l - z) - (1 + 16 z 2 + ^ 4 ) (6 Li(s) - tt 2 
-32 2 3 / 2 (l + z) [vr 2 - 4 Li(y/z) + 4 Li(-^i) - 2 ln(s) In 



(A.l) 



In the counterterms to the two-loop matrix elements we use the following functions: 



B(a) 
B x (a) 

B xx (a) 
B 2 (a) 

B 2 x{a) 

B2xx{o) 



ln(l — x(l — x)a) 
xln(l — x(l — x)a) 
x 2 ln(l — x(l — x)a) 
(ln(l -x(l -x)a)f 



[ x(\n(l-x(l-x)a)Y 
Jo 



x 2 (ln(l -x(l -x)a)f 



(A.2) 



For the explicit expressions one should check a = a + Ie in order to remain on the correct 
Riemann sheet. Further functions are needed for the unexpanded counterterms: 



C(a 
C x (a 

C xx ( a 
C 2 (a 
C 2x {a 



2(x(l - x)a)/(l-x(l-x)a) 

2x(x(l — x)a)/{l — x{l — x)a) 

2x 2 (x(l - x)a)/(l - x(l - x)a) 

41n(l — x(l — x)a)(x(l — x)a)/(l — x(l — x)a) 

4a;ln(l — x(l — x)a)(x(l — x)a)/(l — x(l — x)a) 
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C2xx(a) = / 4x 2 ln(l — x(l — x)a)(x(l — x)a)/(l — x(l — x)a) 
Jo 



(A.3) 



Appendix 2: Effective Wilson coefficients 

The effective coefficients C?_ w appearing in Eq. (J2.7|) are defined in our notation as, 

An „ , , 1„ , x 4„ , , 20 „ , , 80^ . , 



where 



I 4 2D 
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— \n(z) H 1 

9 y J 27 9 s 
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Az 



Az 



2 arctan 



4z-s 



V y/s-Vs-4; 



I TC 



for s < A z , 
for s > 4 z . 



(A.4) 



(A.5) 



Note that specific one- and two-loop and matrix-element contributions of the four-quark operators 
0i_6 (including the corresponding bremsstrahlung contributions) such as the one shown in Fig. IA1I 
are included by the redefinition of the Wilson coefficients CV, C 9 and C\ 9 given in (1A~4|) . In fact, 
using this redefinition, the bremsstrahlung and virtual corrections that are shown in Fig. IA2I (see 
next subsection) automatically take these effects into account. The Wilson coefficients Ci in (jA.4|) . 
which are needed to NNLL precision, are presented in [TT|. IU| . For completeness we quote them here 
again: 



Appendix 3: IR virtual and bremsstrahlung corrections 

The universal 0(a s ) bremsstrahlung and the corresponding infrared (IR) virtual corrections which 
can be absorbed into the Wilson coefficients (see (|2.7j) ) were calculated in [121111120] and given by 

a 9 (s) = a(s) + - , a 7 (s) = a(s) + ~ - | In (— 

2 6 3 Vmft 
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H = 2.5 GeV 


/i = 5 GeV 


H = 10 GeV 


a s 


0.267 


0.215 


0.180 


c[ 0} 


-0.697 


-0.487 


-0.326 




1.046 


1.024 


1.011 


(C C ) 


(-0.360, 0.031) 


(-0.321, 0.019) 


(-0.287, 0.008) 




-0.164 


-0.148 


-0.134 




(4.241, - 0.170) 


(4.129, 0.013) 


(4.131, 0.155) 




(0.115, 0.278) 


(0.374, 0.251) 


(0.576, 0.231) 


(PS 1 *, u?>) 


(0.045, 0.023) 


(0.032, 0.016) 


(0.022, 0.011) 


«\ < j ) 


(0.044, 0.016) 


(0.032, 0.012) 


(0.022, 0.009) 


/^eff(O) Pycif(lh 
l°10 > °10 J 


(-4.372, 0.135) 


(-4.372, 0.135) 


(-4.372, 0.135) 



Table 3: Numerical values of Wilson coefficients of (IA.4|I for three different values of /i; the a s expansion of the 
terms are defined by Cf = cf (0) + cf (1) + ... . 




Figure Al: Examples of virtual and bremsstrahlung contributions of the four-quark operators that are taken 

into account by the redefinition of the Wilson coefficients in Eq. (|A.4|I . 



a(s) = -l L i 2 (s)-^ln(s)ln(l-s) 



V - ln(l 

9 V 



-(l-a)ln(l-s) . (A.6) 



The remaining (finite) non-universal bremsstrahlung are encoded in rate (see 12 .5j) ) and FB 
asymmetry (see I2.6|l ). We note that we have chosen the universal functions <7j in ()A.6j) such that 
the non-universal contributions to the rate, namely T77, T99, r 79 , vanish in the limit s — > 1: 



T-n(s) 
T99 (s) 

T71o(s) 

T-gio(s) 



9(2 + s) 
4 

9(1 + 2s) 

\2 



2(l- g ) 2 ln(l- g ) + 6g(2 2 \ 2 S2) Hs) 

(1 - s) 2 

2d -,r H i-H 33(1 ^ m.) 



11 - 7s- 10s' 
(1-s) 
3(1 -3s 2 ) 



4(1 -) 2 ln(1 _ s) _ 4s(3-2s) 2(5-3.) 



9s 



9(1 



9(1 -a) 



1 



2 3(1 -3s) 3 
5 1 



1 s(6 - 7s) ln(s) l(3 -7s + 4s 2 )ln(l-s) | / 7 (s) 



(1 



9 



2 3(1 -s) 3 ((1 



1 s(6 - 7s) ln(s) 2 (3 - 5s + 2s 2 ) ln(l - s) / 9 (s) 



9 



(A.7) 
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Figure A2: Virtual (up) and real (down) QCD corrections generating the terms r, and tr^ 



-{24(1 + 13s - 4s 2 )Li 2 ( v / s) + 12(1 - 17s + 6s 2 )Li 2 (s) + 6s(6 - 7s) ln(s) 



/9(a) 



6(s - l) 2 

+24(1 - s) 2 ln(s) ln(l - s) + 12(— 13 + 16s - 3s 2 )[ln(l - ^/S) - ln(l - s)] 

+39 - 2vr 2 + 252s - 26vr 2 s + 21s 2 + 8vr 2 s 2 - 180^ - 132sv^} , (A.9) 



{48s(-5 + 2s)Li 2 ( v / s) + 24(-l + 7s - 3s 2 )Li 2 (s) + 6s(-6 + 7s) ln(s) 

(A.10) 



6(s- l) 2 

-24(1 - s) 2 ln(s) ln(l - s) + 24(5 - 7s + 2s 2 )[ln(l - yfi) - ln(l - s)] 
-21 - 156s + 20vr 2 s + 9s 2 - 8tt 2 s 2 + 120^ + 485^} . 



Appendix 4: Complete set of scalar integrals 



In (26| it was shown that there are ten linear combinations of the integrals 

(p ■ k) a (q ■ k) b 



Cti Ct2 «3 



d n pd n q 



[{p + k) 2 + m 2 ] ai (q 2 + m|) Q2 (r 2 + m? 3 



\CX3 



(A.11) 



with a + 6 < 3, which are sufficient for treating a// two-loop Feynman diagrams which one can 
encounter in renormalizable theories. Because the ultravilet behaviour of the functions Tii is loga- 
rithmic only, one finds simple finite integral representations: 
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n 3 



71 



1 1 / \ 1 7T 2 
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(A.12) 



d = 4 — 2e is the space-time dimension, and 7 m = 7 + log(7rm 2 //z 2 ). The special functions hi which 
appear in the formulae above are the finite part in the 1/e expansion of Hi and cannot be further 
integrated into well-studied functions, such as the familiar polylogarithms: 



fox (mi, m,2, m 3', k 2 
h 2 (m 1 ,m 2 ,m 3 ; k 2 
h 3 {m 1 ,m 2 ,m 3 ] k 2 
h^rrii, m 2 , m 3 ; k 2 
h 5 (m 1 ,m 2} m 3 ; k 2 
h 6 (m 1 ,m 2l m 3 ; k 2 
hr(mi, m 2 , m 3 ; k 2 
h 8 (m 1 ,m 2 ,m 3 ; k 2 
h 9 (m 1 ,m 2 ,m 3 ; k 2 
h w (m 1 ,m 2 ,m 3 ; k 2 



dxg(x) 
dx [g(x 
dx [g(x 
dx [g(x 
dx [g(x 
dx [g(x 
dx [g(x 
dx [g(x 



dx \q(x 



dx [g(x 



+ fi(x)} 

+ f 1 (x)](l-x) 

+ fi(x) + f 2 (x)} 

+ h(x) + f 2 (x)](l-x) 

+ }\(x) + f 2 (x)](l-x) 2 

+ Ux) + Ux) + Ux)] 

+ fi(x)+f 2 (x) + f 3 (x)](l-x) 

+ fi(x) + f 2 (x) + f 3 (x)} (1 - xf 

+ Ux) + Ux) + Ux)](l-xf 



(A.13) 



The four building blocks g(x), fi(x), f 2 (x), and f 3 (x) of these one-dimensional integral representa- 
tions are explicitly given in section l3~2l 



Appendix 5: Anomalous dimensions 

We quote here all anomalous dimensions necessary for the calculation. They were presented in 
[T?| HHj and have been checked by us. 
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The counterterm contribution which are proportional to C\ and C 2 due to the mixing of the two 
operators 0\ and 2 with the operators Oj (j = 1...12) are given by the following matrix elements 

12 

(s£ + r\ ]T SZyOj \b), i = l,2, (A. 14) 

i=i 

with the renormalization constants = 5ij + SZij where 

= (< + J<$) + ^ (< + 7< + ^4) + 0(c{). (A.15) 

There are two evanescent operators On and 0\ 2 involved in the mixing in addition to the operator 
basis given in 1)2.3)1 . As usual, their choice is not unique, but we follow the ones used by the authors 
of [T7j so as to be able to use their results for the Wilson coefficients: 

O u = (HWvlvT a c L )(c L i l i , i T T a b L )-\S0 1 , 

12 = {s Ll ^ ula c L ){c Ll » 1 v 1 °b L )-\§0 2 . (A. 16) 

The anomalous dimensions are given then by: 
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Appendix 6: 0(A^ r -n/m?.) corrections. 



L QCD/ " L c 

The explicit expressions of the 0(AQ CD /m^) non-factorizable corrections to R(s) and Aps are 



5 1/m 2R(s) 



8A 2 a 2 em V; s V ts (l-s) 2 
9m 2 c 4-K 2 V; b V tb f(z)<z) 



Re 



1 + 6s - s 2 



4z 



C 2 C™ w (s)* 



newf y 



A 2 a 2 em V*V t s (1-sf 
3m 2 c 4vr2 V; b V tb f{z)K{z) 
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Rc 



:i + 3s)F 



43 



c 2 c™ w ( s y 



(A.17) 
(A.18) 



where 



F(r) 
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rfl — r) 
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(A.19) 
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